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Abstract. 

At large distances and in the low temperature phase, the quenched correlation func- 
tions in the 2d random phase sine-Gordon model have been argued to be of the form : 
( [ip(x) — ^(O)] 2 )* = ^4(log |x|) + £>e 2 (log \x\) 2 , with e = (T — T c ). However, renormalization 
group computations predict B ^ while variational approaches (which are supposed to be 
exact for models with a large number of components) give B = 0. 

We introduce a large N version of the random phase sine-Gordon model. Using non- 
Abelian bosonization and renormalization group techniques, we show that the correlation 
functions of our models have the above form but with a coefficient B suppressed by a factor 
1/N 3 compared to A. 
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The 2d random phase sine-Gordon model has been introduced to describe many dis- 
ordered systems including the 2d XY model in a random field [|J, interfacial roughening 
transition 0, randomly pinned flux lines in supraconductors ||, etc.... Its action is : 

The coupling K is proportional to the inverse temperature, K oc 1/T. In addition to the 
random phase, we also introduced a random potential since this is required by one-loop 
renormalization IE. The quenched random variables A v (x) and £(x) are gaussian with law 



-(c^) 2 ~ A v {x)d v y - £(x)e i<fi - C(x)t 



-up 



P[A] = exp 



1 

2^ 



d 2 x 

4:71 



and P[£] — exp 



1 

2a 



d 2 x 
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There are two different phases : a low temperature phase K > K c = 1 in which the 
disorder is relevant, and a high temperature phase K < K c = 1 in which the disorder is 
irrelevant. At K — K c , the disorder is marginally irrelevant and only induces logarithmic 
corrections to the pure system. Of particular interest is the large distance behavior of the 
quenched average of the correlation functions of <f(x). The first proposal for this correlation 
in the low temperature phase was found using renormalization group (RG) techniques to be 
of the following form |], f|, [5|: 

( - ¥>(0)H, = A(\og \x\) + Be 2 (\og \x\) 2 (3) 
with B 7^ and e = ( K ~k Kc ) ^ 1- There is a crossover from a (log |x| ) to a (log|x|) 2 



^2, for e 



(^) 



< 1. 



behavior. The crossover length R cross is such that (logi? c 
It is exponentially large close to the phase transition. 

However the formula (0) is still controversial, theoretically as well as numerically. In 
fact variational approaches predict B = and a (log |x| ) behavior @, 0? HI- The RG flows 
was also found to be unstable again assymmetric replica perturbations @. The numerical 
verifications of ([|) are also not settled : a (log |x|) behavior was found in ref. [10, 11], while 
more recent simulations |TJ| seem to indicate a (log |x|) 2 behavior. In view of this conflict, 
and since the variational approaches are argued to be exact for systems with a large number 
of components we studied a large iV version of the model ([!]). As explained below, using 
RG computations based on the (a priori symmetric) replica trick, we find that our large iV 
model possesses a non-trivial infrared fixed point in which the correlation functions have 
the form (§) but with a coefficient B suppressed by a factor (1/iV 3 ) compared to A. The 
occurence of this factor could explain why the (log |x|) 2 term does not manifest itself in the 
variational approaches. 

•) The large N bosonic and fermionic actions. 

To introduce the large iV version of ([!]), it is convenient to fermionize it. The fermionic 

form of the random phase sine-Gordon model is a massless Thirring model coupled to a 

quenched potential A^ and a random phase £. To define its large N version, we need to 

k 

introduce iV Dirac fermions with components and ip± with k — 1, • • • , N. Let z = x + iy 

and z = x 



iy be the complex coordinates on the plane. The action is : 



S (N) 



/?(&•- 



■ k d^ k , + if). 
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In absence of disorder, it is conformally invariant. The random potential = (A-, A 

coupled to the U(l) currents J z = X)fcLi and J z = J2k=i '4>--k'4>+ °f the unperturbed 

theory. At £ = 0, the random potential does not break conformal invariance. 

There are a priori many ways to generalize the action ([I]) to a large N version. The action 
(fD has been designed in such way as (i) to keep the number of disordered variables fixed, 
(ii) to preserve the exact conformal invariance in absence of disorder, and (iii) to parallel as 
much as possible standard properties of large N models. In particular, we can implement a 
Hubbard-Stratonovitch transformation in order to disentangle the quartic interaction in (f|). 
The action is then quadratic in the N fermions and we can integrate over them. This leads 
to the following representation of the partition function at fixed disorder : 



J exp (-SW) = J DQexp -N J ^Q Z Q- - Tr (log if 



(5) 



with H 



Therefore, at large N the path integral (§|) is 



£ d z ~~ A z — Q 2 

^-a z -q- e 

dominated by the saddle point as usual in large N technique. However, this is not the way 
we will pursue. 

The fermionic action (^) can be bosonized back using non-Abelian bosonization [|14j] . As 
usual, since the pure system describes N Dirac fermions, the pure bosonized theory will be 
described by a su(N) Wess-Zumino-Witten (WZW) model at level one plus a massless free 
field. The su(N) WZW model at level one possesses primary fields taking values in the 
(N — 1) fundamental representations of su(N). Let <^ and a; fc be the chiral WZW primary 
fields which take values in the defining representation of su(N) and in its complex conjugate. 
Their conformal weights are both equal to (^r) • Let us denote by <p the gaussian free field. 
The original fermions ip± can be written as the product of these WZW primary fields by a 

', and 



vertex operator of the gaussian model. Namely, ip+ 

similarly for the other chiral components t/> + and ip_ 
is : 



and ip-- k = a;fc e 
The bosonic form of the action 



S (N) 



s„ 



+ 



K rd 2 x 
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A v (x)d v <p + £(x)$ + C(x)$* e"^ v 



(6) 



where S wzw refers to the su(N)i WZW action. We have introduced the composite fields 



Hz,z) = ELi^ui^kiz) and $* = ELi^k(z)^(z). Equivalent^, $ = tr D (G) with G 



the group valued field of the WZW model. For N = 1 the WZW terms are absent and we 
recover the action ([!]) of the random phase sine-Gordon model. We used the action (^) as 
definition of the model. 

The effect of the Thirring interaction, specified by the coupling constant a, is summarized 
in the normalization of the gaussian action for ip. That is, if is a known function of a whose 
explicit form is not needed. The coefficient K fixes the scale of the dimension of the vertex 

operators : dim ^e^v 1 ^ = Since the dimension of $ is (^r")> the dimension of the 

field coupled to the random variable £ is : 



From the Harris criterion we know that the disordered variable £ will be relevant if this 
dimension is less than one. Thus, the critical temperature is K c — 1 : with our convention 
the critical temperature is independent of N. For K > K c the disorder is relevant while for 
K < K c it is irrelevant. 

The WZW model possesses chiral conserved currents, which we denote by J%(z) and 
J- a (^z) with a = 1, • • • , dim su(N). In order to fixe the normalization in the WZW sector for 
later convenience, we give the operator product expansions of the primary fields n and (j) g 
and of the currents : 

J xab ■ fabc 

J a Az)rM = °— + J— r M {o) + - (8) 

■£(*)#(0) = ^^(0) + --- (9) 
#(*)M0) = -Lr($ + *tffj h + (10) 

where [t a ,t & ] = —if abc t c with / a6c the structure constants of su(N) and tr{t a t b ) = I n S ab . 
Eq.(pp encodes the fact that the level of the representation of the su(N) affine algebra is 
one. We will also need the Casimir operator whose value in the adjoint representation 
is denoted by Cq and in the defining representation by C n . We have NC n = I n dc with 
do = dim(su(N)) = N 2 — 1. A convenient normalization is Co = 2N, C n = (N 2 — 1)/N and 
thus L = 1. 



Similarly as the random phase sine-Gordon model ||15|| , the large N model possesses a 
U(l) symmetry which amounts to absorbe any translation of <p into the random variables 
and £. Its Noether current is represented by insertions of (d^tp) in the quenched connected 
correlation functions. This symmetry implies that all the quenched connected correlation 
functions of (d^ip) are unaffected by the disorder. But this symmetry has another conse- 
quence : the ^-dependence of the quenched correlation functions can be factorized. The 
simplest way to visualize it consists in noticing that the A^-dependence of the action can 
be absorbed into a shift of if. Indeed, let us decompose as A^ = d^A + e^ u d u (, (this is 
always possible on the plane). The field ( decouples from the action and from the measure. 
It is therefore irrelevant and we can set it to zero. Denoting by S^((p\A; £) the action (j^j 

1 r d 2 x 



with A^ = d^A, we have : 



S^MA, = - A/K\A = 0, Q - ^ J ^(A)) 
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with | = ^e tk l K . For the correlation functions involving the vertex operators e tapip( - Xp \ this 
implies : 

Let G ai> ...(xx, ■ ■ • \g, a) be their quenched averages. Using the fact that £ and £ have the 
same measure, and integrating (p|) over A using the free field gaussian measure (H), P [A^] = 



exp 



~ J ^(d^A) 2 with A^ = d^A, we deduce that : 

G aij ...( Xl , ■ ■ ■ \g, a) = J] K ~ x q \ 29a ^ K2 G ait ...(xi, ■■■\g = 0,a) (12) 
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Equivalent ly, 



d g G aij ...(xi,---\g,cr) = I XJ-^logflxp - x q \ 2 )J G aij ...(xx, ■ ■ ■ \g, a) (13) 

This Ward identity will be useful for analyzing the renormalization group equations. 
•) The effective action and the beta functions. 

To compute the renormalization group equations we used the (a priori symmetric) replica 
trick. Introducing n copies of the system with the same disorder and then averaging over 
the disorder gives the following effective action : 

j («" ** ■ e**-*>) - A / (XVfT j (14) 

with G rs = K5 rs — g. In eq. fll4j) , the indices r, s, ■ ■ ■ refer to the replica index and therefore 
run from 1 to n. The fields with a replica index refer to the copies of the original fields in 
the r th replicated system. The extra A-term arises from the regularization of the cr-term for 
r = s. Even if we would not have introduced it at this point, it would have been generated at 
one-loop. Since it does not couple the replica, we could have introduced it in eq.(||). There 
it represents a current- current interaction in the WZW sector. 

Since the interaction only involves the difference of the fields (p r , it is convenient to 
decompose their kinetic term as follows : 

l -{d^)G rs {d^ s ) = {K ~ n n9) (^(I> r )) +^E0W-^)) 2 (15) 

The field (J2 r ^^ r ) decouples. Since its correlation functions represent the averages of 
the connected correlation functions of the U(l) current, we recover that these correlations 
are unaffected by the disorder. The Ward identity ([T3| ) can also be recovered using this 
decomposition. 

The renormalization group allows us to perturbatively analyse the behavior of the system 
in the low temperature phase. We present a one-loop computation, which is valid close to the 
critical temperature, ie. ( K ^ Kc ) ^ 1- We recall that if the partition function and correlation 
functions are computed with the measure / D(p exp(— S) with S = + J2i g % I d 2 xOi(x), 
where S 1 * is the unperturbed fixed point action and Oi(x) are a set of relevant primary 
operators of dimension hi, the beta functions at one- loop are : 

g* = (3%g) = (2 - hjg* - vr £ C)Jg k + ■■■ (16) 

jk 

The summation in the second term is over all the relevant fields generated by the operator 
product expansions, the coefficients Cj A are determined by the relation : Oi(x)Oj(0) = 
J2 k \x\ h *- hi - h iCtO k {0) + --- 
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Let us introduce the following notation for the perturbing fields in the action (|i4|) : 

Q 1 = E ( $r $ * ' e^-^) (17) 

°* = ^E^C^-^W-v*) ( 18 ) 



03 = EfE^ a;r ) 

r V a / 



(19) 



The fields 0\ and C 3 are the perturbing fields associated to the coupling constants a and 
A, respectively. The field O2 is one of the kinetic field for the </9 r 's. We have the following 
operator product expansions : 

m , ^ / s 2JV(n-2)^ , v 2Nn m . , 2N(n-l)/I U/n , s 

= ^ i(O) + ]^=2^( ) - \ z \2 J ° O 3 (0) + ■ ■ ■ 

03(^3(0) = -^|O 3 (0) + ... (20) 

\ z \ 

OaWOi(O) = --^O 1 (0) + ... 

l z l 

where the dots refer to irrelevant terms. The operator product expansions 02(^)^2(0) and 
Oz(z) ©2(0) only contain irrelevant terms. Using the formula ([16]) we get the beta functions 
at finite n. Since the field (J2r V^) decouples, the coupling (K — ng) is unrenormalized at 
any order in perturbation theory, and (3k = n/3 g . Setting n = as required by the replica 
trick, we get for (^f^) < 1 : 

f3 a = l(^J^) a -2Na 2 -C n Xa + --- 



C Gx2 2N 2 



A = -iT^-T 7 ^" (21) 

and /?x = 0. Here the dots refer to higher loop contributions. So, K is unrenormalized at 
n — 0. These equations show that the couplings to and to O3 are generated at one-loop 
even if one does not include them at tree level. Notice the opposite signs in the A 2 and a 2 
contribution to (3\. 

From equation (|2lD, we immediatly see that in the low temperature phase (K > K c ), the 
beta functions (3 a and f3\ possess a non-trivial infrared zero at cr* and A* given by : 

R Kc ) = 2Na* + C n K~ N ^ 00 4N a* (22) 

The last beta function f} g does not vanish at a*, A*, but it behaves like : 
1 (K-K c \ 2 K-K c 

^-smi^r) + "'- Sor N>>1 ' ~kT < 1 (23) 
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Hence, even at a = ct* and A = A*, the coupling g will continue to flow. We may characterize 
such pseudo-fixed point as a "run away fixed point" . 



•) The RG equations and the quenched correlation functions. 

The non-vanishing of the beta function f} g at the infrared pseudo-fixed point has direct 
consequences. Consider the quenched correlation functions G ai) ...(ii, • • •) involving the vertex 
operators e iapip ^ Xp \ They satisfy the renormalization group equations : 



d d 

Efc + E7 P - E /w 

P Udj p P c=g,cr,X UU 



G au ...(xi, ■ ■ ■) = 



Here 7 P are the anomalous dimensions. Since the (^-dependence is explicitly known from 
eq. (JITf) , we can separate the function j3 g from (3 a and (3\ to get : 



v u ^p p c=<t,a uu p<q 



G au ...(xi, ■ ■ •) = 



In particular at the infrared fixed point cr*, A*, in which 3* = /3{ = 0, we get : 



P pt T v 
p ^^p p 



p<q 



Gl,...(xi,...) = 



(24) 



where 7* are the values of the anomalous dimensions at the infrared fixed point and f3* is 
given in eq.(p3|) for large N and for ( K x Kc j "C 1. The net effect of the g-flow is to add the 
extra logarithmic term in the renormalization group equations (|24]) . 

Eq.(^) can be used to compute two-point functions at the infrared fixed point. Con- 
sider : 



Gi(x) = (exp (ia((p(x) -<p(0))) t 
G 2 (x) - 



= (exp (iai(p(x))) (exp (-ia¥>(0)))„ 

Let 7* 2 be their anomalous dimensions at the infrared fixed point. The renormalization 
group equation (E3) gives : 



G 1>2 {x) = \x\ 



exp 



2K 2 



(log|x|) ; 



(25) 



Notice the (log|x|) 2 correction which arises from the renormalization of g. The anomalous 
dimensions are even in the charge a and vanishes at a = 0, therefore : 



7l,2 = ^Pl,2 + £V), With P \ = l + 



( K-K c 



(26) 



Expanding eq.(KH) in power of a 2 gives the two-point functions of <p 



[ (<P(x) - <P(0)) \\ 



2p\ 


log 


\x\ + 


P9 

2K 2 


K 


2p*2 


log 


\x\ + 


3* 
2K 2 


K 



(27) 
(28) 



Comparing the formula for p\ and (3* close to the phase transition, eq. (^6|,^|), we 
see that at large N the (log|x|) 2 term is suppressed by a factor of 1/iV 3 compared to the 
(log|x|) term. Hence there is no a priori contradiction between the replica symmetric RG 
computations and the variational approaches : they coincide in the infinite N limit, i.e. in 
the regime where the latter are expected to be exact. It would be interesting to compute 
the -k corrections in the variational approaches for comparison. 

Note that the (log \x\) 2 cancel in the connected correlation function ([<p(x) — (p(0)] 2 ) conn as 
it should be, since this connected correlation function is unaffected by the disorder. Eq. 



28|) are exact to all orders in perturbation theory provided we may trust the renormalization 
group using symmetric replica. Only the exact value of the beta function f3* depends on the 
order of the perturbation theory. 

Acknowledgement : It is a pleasure to thank Henri Orland for discussions. 

References 

[I] J.L. Cardy and S. Ostlund, Phys. Rev. B 25 (1982) 6899. 
[2] D. Cule and Y. Shapir, Phys. Rev. Lett. 74 (1995) 114. 
[3] M.P.A. Fisher, Phys. Rev. Lett. 62 (1989) 1415. 

[4] J. Villain and J.F. Fernandez, Z. Phys. B54 (1984) 139. 

[5] J. Toner and D. DiVincenzo, Phys. Rev. B 41 (1990) 632. 

[6] S.E. Korshunov, Phys. Rev. B 48 (1993) 3969. 

[7] T. Giamarchi and P. Le Doussal, Phys. Rev. Lett. 72 (1994) 1530. 

[8] H. Orland and Y. Shapir, Europhys. Lett. 30 (1995) 103. 

[9] T. Giamarchi and P. Le Doussal, Phys. Rev. Lett. 74 (1995) 606. 

[10] G. Batrouni and T. Hwa, Phys. Rev. Lett. 72 (1994) 4133. 

[II] H. Rieger, |cond-mat/9502092| . 



[12] E. Marinari, R. Monasson and J. Ruiz-Lorenzo, |cond-mat / 9503074 



[13] M. Mezard and G. Parisi, J. Phys. A 23 (1990) L1229. 
[14] E. Witten, Comm. Math. Phys. 92 (1984) 455. 

[15] U. Schulz, J. Villain, E. Brezin and H. Orland, J. Stat. Phys. 51 (1988) 1. 



8 



